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We disclose the origin of anisotropic percolation perimeters in terms of the Stochastic Loewner 
Evolution (SLE) process. Precisely, our results from extensive numerical simulations indicate that 
the perimeters of multi-layered and directed percolation clusters at criticality are the scaling limits 
of the Loewner evolution of an anomalous Brownian motion, being subdiffusive and superdiffusive, 
respectively. The connection between anomalous diffusion and fractal anisotropy is further tested by 
using long-range power-law correlated time series (fractional Brownian motion) as driving functions 
in the evolution process. The fact that the resulting traces are distinctively anisotropic corrobo¬ 
rates our hypothesis. Under the conceptual framework of SLE, our study therefore reveals new 
perspectives for mathematical and physical interpretations of non-Markovian processes in terms of 
anisotropic paths at criticality and vice-versa. 
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The Stochastic Loewner Evolution (SLE) [l| has rev¬ 
olutionized our understanding of two dimensional loop¬ 
less paths, as recognized among others by several Fields 
medals UM- It provides a mapping between these paths 
and a real valued function, called “driving function”, that 
is a random walk if the path is a conformally invari¬ 
ant fractal. This establishes a relation between the frac¬ 
tal dimension of the path and the diffusion constant of 
the random walk. Although several generalizations have 
been proposed , due to its nature, SLE has been re¬ 
stricted to isotropic models dSl. However, anisotropic 
paths, namely, paths with a preferential direction, appear 
quite commonly in Physics. By numerically determin¬ 
ing the driving function of anisotropic paths, we discover 
that they are consistently mapped onto correlated ran¬ 
dom walks, meaning that the Markovian property of the 
driving function is violated. More precisely, we show that 
resulting anomalous diffusion is characterized by an ex¬ 
ponent that is related to the degree of anisotropy. This 
behavior can be subdiffusive, as it is the case for the 
hull of directed percolation or superdiffusive, as found 
for multi-layered percolation. 

The chordal variety of SLE (the one we will focus on 
this work) deals with curves y* that start at the origin and 
grow towards infinity while restricting themselves to the 
complex upper half-plane H. The curve yt is connected 
to a real-valued driving function Ut through the relation 
yt = Ot^iUt) where gt{z) is the solution of Loewner’s 
equation 3 23| . 


dtgt{z) = - -—, go{z)=z. (1) 

gt[z) - Ut 

In his seminal work [ij, Schramm showed that if the 
measure over yt displays conformal invariance and do¬ 
main Markov property, then the only possibility is that 
Ut be a Brownian motion with a single free parame¬ 


ter K, the diffusion coefficient. This is often written as 
Ut = y/KBt where Bt is a standard Brownian motion 
(with diffusion coefficient equal to unity). The value of 
K is related to the geometric properties of yt, includ¬ 
ing its fractal dimension, which is determined by the re¬ 
lation df = min(l-|-|,2) [^. A few lattice models 
have been shown to converge to SLE in the continuum 
limit ll. Il2l l22l. and many more are conjectured to do 
so [^, ll.ll. Il4 . Of particular interest is the proof 
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FIG. 1. Percolation models used to generate the SLE curves, 
(a) Regular percolation, where each site is occupied with the 
same probability p [l^. (b) Multi-layered percolation, where 
some rows are occupied with probability p + A (dark gray 
rows) and others with p — A (light gray rows) [1(| . (c) Di¬ 
rected percolation is a spreading process which starts at the 
bottom of the tilted lattice and can only advance upwards 
with probability p 
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FIG. 2. Examples of cluster perimeters at the critical point of (a) isotropic percolation on a triangular lattice (pc = 0.5), (d) 
multi-layered percolation also on a triangular lattice (A = 0.2 and Pc = 0.5), and (g) directed percolation on a square lattice 
(pc ~ 0.644). A detail of each curve in (a), (d) and (g) (pink square) can be seen in (b), (e) and (h), respectively. The driving 
fnnctions obtained by applying the zipper algorithm to the curves (a), (d) and (g) are shown in (c), (f) and (i), respectively. 


that the perimeter of a percolation cluster on a triangular 
lattice follows SLE with /t = 6 as a scaling limit, allowing 
for a formal computation of its critical exponents Q. 

In this work we explore the possibility of using Loewner 
evolutions to study anisotropic fractal systems, i.e., sys¬ 
tems with different critical exponents in each direction. 
These systems are not scale invariant, therefore not con¬ 
formally invariant either. We are particularly inter¬ 
ested in two variants of the percolation model that show 
anisotropic behavior, namely, multi-layered percolation 
and directed percolation (see Fig. [T|). Precisely, we gen¬ 
erate the border of percolating clusters, numerically com¬ 
pute their corresponding driving function, and then ana¬ 
lyze the diffusive properties of these numerical sequences. 
In the general case, we expect that the mean square dis¬ 
placement of Ut behaves like, 

(u?) ^ br ( 2 ) 

as t —>■ oo. In the case of traditional SLE, a = 1 and 
b = K. The driving functions of anisotropic percolation 
models, we found, display very distinctive anomalous dif¬ 
fusive behavior (a ^ 1). Finally, we show that our ap¬ 
proach is also valid in the opposite direction, namely, 
the SLE consistently leads anomalously diffusive driving 
functions to traces that display clear anisotropic scaling. 

In order to evaluate the SLE driving function of the 
cluster perimeters we used the zipper algorithm with a 
vertical slit discretization [3) [2^. In this method, given 
a lattice curve { 0 , 71 , 72 ,..., tat}, its driving function can 
be recovered by applying the relations, 

^ fc k 

4 = - ^ Im{wj }2 Ut^ = Re{a;j}, (3) 
i=i i=i 


where the cokS are determined recursively by 

Wfe =/fe-i o/fc _2 o ... o/i( 7 fc) a;i= 7 i, (4) 

and fkiz) 

fk{z) = i\J-lv[i{LOkY - (z: - Re{uJk}f. ( 5 ) 

This algorithm, however, does not guarantee that the dis¬ 
cretized times tk are equally distributed, even for curves 
of same length and step size. To obtain an ensemble of 
curves defined for the same time sequence, we linearly in¬ 
terpolate the obtained driving function at equally spaced 
points in logarithmic time in the interval [l,logf/], for 
some suitable tf. 

As already mentioned, we also performed the oppo¬ 
site operation of computing the SLE trace from a given 
driving function The process is simply the inversion of 
the algorithm previously described. Given a discretized 
driving function Ut = { 0 , Ut -^,..., t/t„}, the trace can be 
obtained by repeatedly applying the functions 

7i = 5o offi o • ■ • offi(O), (6) 

where the mappings are also chosen to represent a vertical 
slit discretization, 

gi{z) = i\j 4(tj - ti-if - + (Ut, - Uti^J. (7) 

Instead of using an approximate algorithm [^, we 
chose to use the one described above, as they are ex¬ 
act (for a given discretization). Their complexity scales 
as 0{N‘^) which can get quite time consuming, specially 
for large values of k, requiring a large number of points 
to get accurate results. We resorted to GPU paralleliza¬ 
tion (where each 7 ^ is computed by a single thread) to 
achieve satisfactory accuracy. 
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FIG. 3. Mean squared displacement of the driving functions 
for the three percolation models studied. The curves are the 
results of the numerical procedure described in the text ap¬ 
plied to 10^ realizations of each type of percolation model. 
The 95% confidence intervals were bootstrapped over 400 re¬ 
samplings [^ . but, being smaller than the symbols, are not 
shown. As expected, in the case of isotropic percolation, the 
displacement scales linearly with time, while it shows instead 
a distinctive subdiffusive behavior for directed percolation, 
with an exponent a « 0.67. In the case of multi-layered per¬ 
colation, a clear superdiffusive behavior, with an exponent 
a « 1.78, can be observed for A = 0.4. The inset shows how 
this anomalous diffusion regime is gradually achieved as we 
increase the degree of anisotropy A. 


We start by testing our approach on standard isotropic 
percolation, which has been extensively studied as a sim¬ 
ple but rather rich and illustrative model for critical¬ 
ity [l^. It is basically a lattice model with binary disor¬ 
der, where each site (or bond) is occupied with probabil¬ 
ity p. For a given critical probability pc, the presence of a 
giant spanning cluster is detected. In the thermodynamic 
limit, if p < Pc, where Pc is the percolation threshold, the 
system never percolates, otherwise it always does. In 
particular, m = 1/2 for site percolation on the triangu¬ 
lar lattice [15|. It has been mathematically proven that 
the perimeter of the giant cluster at the critical point 
follows SLE with k = 6 Q . We perform simulations with 
10 “^ realizations of percolation perimeters of length 10® 
lattice units generated using the algorithm described in 
Ref. [s^ on the triangular lattice. Fixed boundary con¬ 
ditions are adopted, in which every site on the left side 
of the bottom row is always unoccupied and the ones on 
the right side are always occupied. In Fig. [5^1 we show 
a typical realization of an isotropic percolation perime¬ 
ter and the corresponding driving function, as computed 
using the algorithm Eqs. [3][5l Finally, from the driving 
functions, we calculate their mean squared displacement 
(t//) as a function of time, and find that k = 6.27 ± 0.30 
and a = 0.996 ± 0.005, as shown in Fig. |3^. 


Unlike regular percolation, where each site or bond is 
occupied with probability p, in multi-layered percolation 
this is done with probability p± A, where A G [0, ^], and 
the signs plus or minus are chosen randomly with equal 
probabilities for each row of the lattice [1^ [Slj. Here, 
the parameter A represents the degree of anisotropy of 
the system, with A = 0 being equivalent to isotropic 
(regular) percolation. 

We generate an ensemble of 10'’’ multi-layered perco¬ 
lation perimeters of length 10® lattice units on a trian¬ 
gular lattice at the critical point for different values of 
A. For every value of A we found Pc = 0.5 using the 
cluster perimeter method [s^. As for the standard per¬ 
colation case, after calculating the driving functions, an 
example of which is shown in Fig. [2]3, we compute the 
corresponding mean square-displacement to find that it 
exhibits characteristic superdiffusive behavior for every 
value of A > 0. As can be observed in the inset of Fig. 3, 
however, a long transient behavior is present for small 
values of A before a distinctive power-law behavior is es¬ 
tablished. For A = 0.4, after a short transient, the least- 
squares ht to the data gives a power-law, (U/) = 6t“, 
with b = 10.38 ± 0.68 and a = 1.78 ± 0.01, which extends 
over more than three orders of magnitude. 

Next we investigate the diffusive behavior of driving 
functions generated from directed percolation perimeters. 
As dehned, directed percolation is a spreading process 
where a cluster can only grow along preselected directions 
in a lattice, and each site is occupied with probability 
P [i3- Shown in Fig. [TJ: is a typical realization of a di¬ 
rected percolation perimeter generated on a tilted square 
lattice at the critical point, pc = 0.644700185(5) [s^. Us¬ 
ing this simulation setup, the perimeters of the spanning 
clusters are obtained here using a simple walker algo¬ 
rithm, as illustrated by the red curve in the example 
shown in Fig. Ic. From the ensemble of the generated 
driving functions, once more the resulting mean square 
displacement displays a characteristic anomalous behav¬ 
ior. Precisely, the least-squares fit to the data in the 
scaling region yields subdiffusive diffusion, as shown in 
Fig. 3, with a pre-factor b = 3.74 ±0.07 and an exponent 
a = 0.676 ±0.001. 


The results obtained by the previous analysis suggest 
that the presence of long-range correlations in the driv¬ 
ing function should lead, through the Loewner evolution 
process, to anisotropic fractal traces, and vice-versa. In 
order to test this hypothesis, we analyze the behavior of 
traces driven by stochastic processes exhibiting anoma¬ 
lous diffusion. We choose to use fractional Brownian time 
series generated according to a given Hurst exponent H, 
which is related to the diffusion exponent by a = 2H 341. 

We generated the drive Ut as a fractional Brownian 
Motion with Hurst exponent H and diffusion constant 
b in N time steps U uniformly spaced in the interval 
[0,t/]. In order to simulate fractional Brownian motions 
with reasonable control over the diffusive constant 6, the 
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Davies-Harte algorithm was used [35| . The 7 t^ were com¬ 
puted from Uti using Eq. |6l We then interpolated the 
trace 7 (f) (the same Ja as before, but parametrized by its 
length instead of the Loewner time) in M equally spaced 
points £i G [Ojimax]- This interpolation step was neces¬ 
sary because the zipper algorithm generates discretized 
traces with highly non-uniform step sizes | 7 i — 7 i-i|. Al¬ 
though this does not diminish the intrinsic error of the al¬ 
gorithm, it makes the analysis easier to perform. In order 
to study whether the scaling is isotropic or anisotropic, 
the root mean squared estimation of the displacement of 
the trace was computed in each direction, that is. 


FxiiAe) 



M-i 

Y, [X{£j+i) - 


3=0 


( 8 ) 


where A(£) = Re{ 7 (£)}. Analogously, Fy(iA^) is defined 
taking instead Y{£) = Im{ 7 (£)}. 

Our numerical scheme was applied to three sets of 
times series, each with 100 realizations generated to re¬ 
produce the corresponding properties (in terms of H, 
b and tf) of the driving functions originated from the 
isotropic and anisotropic percolation traces previously in¬ 
vestigated. More precisely, the first set of time series cor¬ 
responds to uncorrelated Brownian motion, the second to 
correlated or persistent, and the third is anti-correlated 
or anti-persistent [s^. The remaining parameters {N, 
M and imax) are chosen to ensure the accuracy of our 
results. The precise values of all parameters adopted in 
the simulations are reported in Table I. Figure |4] shows 
that the traces evolved from these time series have sim¬ 
ilar behavior to their corresponding percolation models, 
in the sense that a clear anisotropy can be observed in 
the correlated and anti-correlated simulations, while the 
uncorrelated one displays isotropic behavior, as expected. 

In summary, our numerical analysis offers compelling 
evidence that a variation of the Stochastic Loewner Evo¬ 
lution, obtained by taking as driving function a stochas¬ 
tic process with anomalous diffusion, may be the scal¬ 
ing limit of anisotropic critical models. In particular 
we looked at two anisotropic variants of percolation: di¬ 
rected percolation and multi-layered percolation. The 
former was found to be associated with subdiffusive driv¬ 
ing functions, while the latter are superdiffusive. We also 
tested the inverse relation, finding that driving functions 
with anomalous diffusion do indeed generate traces with 
anisotropic features. This possibility opens new ques¬ 
tions, like how the critical exponents of SLE traces de¬ 
pend on the addition of long-term correlations to the 
driving function. Moreover it would be interesting to 
know if one can obtain exponents of actual physical mod¬ 
els with such a generalized theory. We expect that the 
further developments of this new variant of SLE may pro¬ 
vide some insight on the critical behavior of anisotropic 
systems, the same way the original SLE was to isotropic 
systems. 



H 

b 

tf 

N 

M 

^max 

Ensemble 1 

0.5 

6.0 

2 ■ 10® 

10® 

10® 

2 ■ 10 ^ 

Ensemble 2 

0.8 

16.0 

3- 10^ 

10® 

10® 

8- 10'‘ 

Ensemble 3 

0.33 

3.8 

5- lO’’ 

10® 

10® 

2 • lO"' 


TABLE I. Simulation parameters used to generate the SLE 
traces. H is the Hurst exponent and b is the diffusion coeffi¬ 
cient of the fractional Brownian motion used as driving func¬ 
tion. The curves were computed for N times ti equally spaced 
in the interval [0, tf]. The resulting trace is reparametrized as 
a function of its length and interpolated in M points equally 
spaced in the interval [0,£m,ax]- 
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